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"in : 

Abstract 

O ; 

The paper studies non-Lie symmetry of the Klein-Gordon-Fock equation (KGF) 
in (p + g)-dimensional Minkowsky space. Full set of symmetry operators for the n- 

■ order KGF equation was explicitly calculated for arbitrary n < oo and p + q < 4. 

■ Definition was given for generalized Killing tensors of rank j and order s, and 
for generalized conformal Killing tensors of rank j and order s as a complete set 
of linearly independent solutions of some overdetermined systems of PDE. These 

■ tensors were found in explicit form for arbitrary fixed j and s in Minkowsky space 
of dimension p + q < 4. The received results can be used in investigation of higher 
symmetries of a wide class of systems of partial differential equations. 

O ' 

in . 

Introduction 

^ : 

^p" 1 ! Classical group theoretical analysis of differential equations whose foundations were laid by 
Sophus Lie over a hundred years ago finds increasing utilization in modern mathematical 
physics (see e.g. [H 121 IS])- At the same time certain limits of the classical Lie approach 
become obvious that nevertheless do not allow full description of the symmetry of an 

> . equation under study [U Ej. In particular, does not allow calculation of higher order 
symmetry operators that are widely used for calculation of reference frames admitting 
solution of equations in separated variables 00 IB], in calculation of motion constants |Hj 
and in many other problems. 

The present paper deals with investigation of non-Lie symmetry of the Klein-Gordon- 
Fock equation in (p + (/)-dimensional Minkowsky space 



where x is a real parameter, 
0, n ^ u, 

<r = { i, n = v< P , (2) 

-1, p<fi = v<p + q, 

<p = (p(xi,x 2 , . . . , Xp+q) is a function of p + q variables. 

1 The first version of the paper was published as a preprint in Russian: Preprint N 90.23, Kyiv, Institute 
of Mathematics, 1990, 59 p. 
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A symmetry operator of equation (1) is an arbitrary operator Q (linear, nonlinear, 
differential, integral) that transforms solutions of this equation into solutions, that is |Hj 

L{Q(p) = 0, if Lip = (3) 

(see Section 1 below for more rigorous definition). 

We will call a differential operator of a finite order n being a symmetry operator of 
equation (1) a n-th order symmetry operator. 

Description of the maximal (in the sense of Lie) symmetry of equation (1) may be 
reduced to finding of all linearly independent first order symmetry operators. Such op- 
erators are well-known, they form a basis of the Lie algebra of the generalized Poincare 
group P(p, q) (for x ^ 0) or for the conformal group in (p + (^-dimensional space (when 
x = 0). 

One of the main results of the present paper is calculation in explicit form of a complete 
set of n-th order symmetry operators of equation (1) for arbitrary n < oo and p + q < 4. 

It is well-known that description of first order symmetry operators is based upon cal- 
culation of explicit form of the Killing vector ^UJ |2| that corresponds to the space of 
independent variables. We associate with higher order symmetry operators more complex 
fundamental objects that we call Killing tensors of rank j and order s, with j, s = 1, 2, . . . 
and conformal Killing tensors of rank j and order s. 

In this paper we give the definition of the mentioned tensors as a a complete set of 
linearly independent solutions of some overdetermined systems of PDE and find these 
tensors were found in explicit form for arbitrary fixed j and s in Minkowsky space of 
dimension p + q < 4. The results can be used in investigation of higher symmetries of 
a wide class of systems of partial differential equations of mathematical physics given in 
the same space, in particular, of relativistic and galilei-invariant wave equations. 

Let us describe briefly arrangement of our presentation. Main definitions related to 
higher order symmetry operators, are adduced in Section 1, definition of Killing tensors 
of rank j and order s is given in Section 2, first order Killing tensors of rank j and order 
s in explicit form are found in Sections 3 and 4, conformal Killing tensors and Killing 
tensors of arbitrary rank and order are shown in Sections 6, 8 and 9. A complete set 
of n-th order symmetry operators for equation (1) with zero and non-zero "mass" x are 
adduced in Sections 5 and 7. 

1 Symmetry operators of order n 

For the purpose of our study it is sufficient to consider only solutions of equation (1) 
defined on an open set D of the four- dimensional manifold M p + g consisting of points with 
co-ordinates (xi, £2, • ■ ■ , x p+q ) and analytical with respect to real variables xi, x 2 , . . . , x p+q . 
The set of all such solutions forms a complex vector space that we designate by the symbol 
Tq. Setting D as fixed (e.g. assuming that D coincides with M p +q), we will call JF the 
set of solutions of equation (1). 

Let us designate with T a vector space of all complex- valued functions defined on D 
and being real analytical, and with L a linear differential operator (1) defined on T . Then 
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Lip G T when ip G T . At that JF is such subspace of the vector space T that coincides 
with zero-space (kernel) of the operator L. 

Let Tl n be a set (class) of differential operators of the order n defined on T . Then a 
symmetry operator Q G 9JT„ of equation (1) is defined as follows. 



Definition. A linear differential operator of order n 

Q = Y,Qi, Qi = H^ a *-^—— 
L — ' ox m 

i=0 



&X ai ^X a2 • • • &'Eai 



(4) 



is called a symmetry operator of equation (1) in the class 9Jl n (or symmetry operator of 
order n) if 



[Q, L] = a Q L, a Q G Wl n - 



(5) 



where [Q,L] = QL — LQ is the commutator of the operators Q and L. 

Relation (5) should be understood in the sense that operators in the right-hand and 
left-hand parts give the same acting on an arbitrary function ip G T . Functions H ai ' a ' 1 ' , ' ,)at 
to be determined are symmetric tensors of rank i. Hereinafter the parentheses enclose the 
set of symmetric indices. 

It is easy to see that the relation (3) follows from (3) for each ip G J-q. The reverse 
statement is also true: if an operator (4) satisfies the relation (3) for arbitrary ifj G Tq then 
the condition (5) is satisfied for such operator with some operator «q. In the case n = 1 
the symmetry operators defined above may be interpreted as generators of the symmetry 
group of the equation being considered jjj. We can show that the set of symmetry 
operators Q G 97ti generates a Lie algebra, and corresponding finite transformations from 
the invariance group may be obtained by integration of the Lie equations (HE]. 

Symmetry operators of order n > 1 are not generators of a Lie algebra anymore and 
characterize generalized (non-Lie) symmetry of an equation under study. The problem of 
description of a complete set of n-th order for equation (1) be reduced to finding of the 
general solution of the operator equations (5). 



2 Equations for coefficients of symmetry operators. 
Killing tensors of rank j and order s 

For simplification of further calculations it is more convenient to present the operator Q 
(4) as the sum of i-multiple anticommutators 



<3 = E4> 

3=0 



(6) 



where 



Qi 



p{a,iai—aj) 



d 

9x ai 







dx 



"2 . 



d 

dx a . 



(7) 
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[A, B] + = AB + BA, F^- a ^ is a symmetric tensor of rank i. Expanding anticommu- 
tators and transferring differentiation operators to the righthand side, it is possible to 
reduce the expression (6) for the operator Q to the form (4), and, vice versa, to write 
down any operator of the form (4) as (6). 

We can use a similar representation for the operator «q e 9Jl n -i 



n-l 



^2 



as = 



=n-2 



/ 



aiO2---0-j 







' dx. 



a i 







' dx 



a 2 







dx n 



(8) 



and write the product cxqL as 



_d_ 

dx. 



+ 



d 
dx» 



1 

+ 2 



a Q 



_d_ 

dx. 



(9) 



where a Qlt = gj. 

Using the representations (6)-(9) and taking into account that 



[Q,L] 



Q, 



d 
dx» 



dQ 



'' dx^ ' 



(10) 



it is possible to reduce the operator equation (5) to the system of equations for coefficients 
ja ia2 ...ai an j pa 1 a 2 ...a j ^ j n f ac ^ substituting (6)-(10) into (5) and putting equal coefficients 
at identical degrees of operators of differentiation, we obtain 



Qfaj+i pa 1 a 2 -a j ) _ 
j(aia2--ai) _ g 



(11) 

(12) 



Here d aj+1 = ^— — , the round brackets contain symmetric indices (so symmetrization is 
implied (11)): 



3 + 1 



_|_ Qa 2 pa 1 a j+1 ...a j _|_ ... _|_ QQj pa\a 2 ...a j+ 1 ^ 



paia 2 ...aj j g a S y mme tric tensor of the rank j. 

If k = 0, the equation for coefficients of the symmetry operator takes the following 
form: 



Qfaj+i pa!a 2 ...aj) _ fi(aja j+1 ja 1 a 2 ...a j - 1 ) 



(13) 



where F a ^- a i and f^^-aj-i are symmetric tensors with zero trace. 

Convoluting equations (13) with respect to one pair of indices, we can eliminate the 
unknown functions j a ^- a ]-^ . As a result we get 



g{aj+l pa 1 a 2 ...a j ) J Qb pb(a 2 a 3 ...a j ga 1 a i+1 ) _ g^ 

m + j — 1 



(14) 



jaia2— Oj-i _ 3 Qb jpb(aia,2—aj-i aj-aj+i) (15) 

m + j — 1 

where m = p + q is dimension of the space of independent variables. 

We see that the problem of description of symmetry operators of order n for the 
equation (1) with x = appears to be equivalent to finding of the general solution of 
the system of partial differential equations given by the formula (11). This system is 
split with respect to the index j, as it splits into independent subsystems corresponding 
to j = 0, 1, ... ,n. As it will be shown below, for complete description of the symmetry 
operators it is actually sufficient to solve only two such subsystems corresponding to j = n 
and j = n — 1. 

In the case j — 1 the system (11) coincides with the Killing equations jSGDj, and for 
j = 2 it coincides with equations for the Killing tensor [Hj in the flat de Sitter space. The 
corresponding equations (14) determine conformal Killing vector and Killing tensor in the 
p + g-dimensional Minkowsky space. 

We shall call functions F a ^- a 3 satisfying equations (11) (or (14)) Killing tensors (or 
conformal Killing tensors) or rank j and order 1. The meaning of the term "order 1" 
(that we will omit sometimes) will be explained below. 

The equations (11), (13) for a tensor of arbitrary rank were introduced (in the case 
j > 2, without relation to any particular problem) in the paper ^Tj. However, the general 
solution of these equations, as far as we are aware, was obtained in an explicit form only 
for j = 1 and j = 2 [T2]. 

In the process of investigation of higher order symmetry operators admitted by systems 
of partial differential equations, we have to deal with more complicated equations for 
coefficients of such operators than those given by formulae (11) or (14). These equations 
include derivatives of the order s > 1 and have the form 



where F a ^- a i is a symmetric tensor, and 



Qa j+1 Qa j+2 . . . Qdj+s pa ia2 . 



SL 



0, (17) 



where F a ^--- a i is a symmetric tensor with zero trace, and the symbol [■] designates 
the zero trace part of the tensor inside the square brackets (in our case it is a symmetric 
tensor of the rank R — j + s) : 

m ( d 

^aia 2 ...a ii j5L _ Q ai a 2 ...a R _|_ ^ ^ ( — j) rf _K" rf j J^J g a 2i-\a.2i 

d=l \i=l 
^ jpa, 2 d+ia2d+2...a i b 1 b 2 b a b 4 ...b 2d _ 1 



9b 1 b 2 9b :i b 4 ...gb 2d _ 1 b 2 



where } is the integer part of the number 



2 > 



Kd (n-2d)\2 d - l ^2(n-l) + m-2 , ^ 
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In the case s = 1 the equations (16) and (17) can be reduced to equations (11) and 
(14) respectively. 

We will call a symmetric tensor F a ^ a z-- a 3 satisfying equations (16) a Killing tensor of 
rank j and order s.We will call a symmetric tensor p a ^- a i with zero trace satisfying 
equations (17) a conformal Killing tensor of rank j and order s. 

In Sections 3-7 below we obtain the general solution of equations (11), (14) for arbitrary 
j in the space of dimension p + q < 4. Equations (16), (17) are discussed in Sections 8, 9 
where their general solution is found for p + q < 4 and arbitrary j and s. 

3 Reduction of equations for symmetry operators 
to a system of linear algebraic equations 

Let us start investigation of the system of equations (11) describing the Killing tensor of 
rank j and order 1. 

The system (11) may be written in the following symbolic form: 

< p»io a ...Oj+i = 0) (20) 

where j^ a ^- a i+^ i s a symmetric tensor of rank j + 1 in m = p + g-dimensional space, and 
unknown functions are components of symmetric tensor of rank j in m-dimensional space. 
Whence we can see that the system under investigation is overdetermined, including ( J ^) 

equations for unknowns, Q) = a \j^i a y designating binomial coefficients. 

Following the general method for solving of overdetermined systems of partial differ- 
ential equations [Tl) . we consider the set of differential consequences of the system (11), 
obtained by differentiation of each term, k times by x^ (i = 1,2, . . . , k). For each fixed k 
such differential consequences are systems of linear homogeneous algebraic equations for 
derivatives 

gbiQb2 _ _ _ Qb k Qa j+1 piaia,2...aj _ p(a\a,2...aj ,dj +1 )b\b2...b k (21) 

These systems have the form 

p(a 1 a2...aj,aj +1 )b 1 b 2 ...b k _ g (22) 

The system of equations (22) determines condition for vanishing of the tensor of rank 
j+k+1 symmetric with respect to j+1 indices a\, 02, . . . , a J+ i and with respect to k indices 
b%, 62, • • • , bk, with unknown components of the tensor (21) of rank j + k + 1 symmetric 
with respect to j indices a\,a 2 , . . . ,a,j and with respect to k + 1 indices Oj+i, 
Whence we conclude that the corresponding numbers of equations (N e ) and of unknown 
variables (iV u ) are given by the formulae 

f j + m\ fk + m-l\ fj + m-l\fk + m\ 

e ~\j+i){ k y N »-{ j )\k+i)> (23) 

where m = p+q is dimension of the Minkowsky space where equations (11) are determined 
(that is the number of independent variables xi,x%, . . . of the function F aia2 - a i) . 
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According to (23) 

N c <N n , k<j, N e = N n , k = j. (24) 

The formulae (23) allow calculation of the number of linearly independent solutions of 
equations (22), as the following statement is true: 

Theorem 1. The system of linear algebraic equations (22) is not degenerate. 
Proof of Theorem 1 is adduced below in Appendix. 

We conclude from (24) in virtue of Theorem 1 that for k = j the system of homogeneous 
linear algebraic equations (22) has only trivial solutions, 

pa\a2...aj,ajj r \b\bi---bj — g 

Whence coefficients of the symmetry operator ~F a ^ a ?-- a i are polynomials on x a (a = 
1, 2, . . . , m) of order j. It follows from (23) that such polynomial contains NJ 1 arbitrary 
parameters, where 




We see that equations (11) have N™ 1 linearly independent solutions that form a com- 
plete system. To find these solutions in explicit form it is necessary to find the general 
solution of the system of linear homogeneous equations (22) for arbitrary given j, m and 
k < j, and then reconstruct polynomials F a ^- a i by found values of derivatives of the 
tensors F a i--- a j> a j+i b i b 2--- b k (i e t us remind that indices after the comma designate derivatives 
with respect to the corresponding arguments). The general solution of equations (11) is 
adduced in Section 4 below. 



4 Explicit form of Killing tensor of rank j 

According to the above proof, calculation of the explicit form of Killing tensor of rank j 
is reduced to finding of the general solution of non-degenerate system of linear homoge- 
neous algebraic equations given by the formula (22). Actual solution of this system with 
arbitrary given j and m is a rather difficult task that may be circumvented using the 
following observation. 

Lemma 1. Let p a ^--- a n-, ft e an arbitrary solution of the system (11) for j = jo, and F a 
be a solution of the same system for j = 1. Then the function 

pa 1 a 2 ...a jQ + 1 _ p(a 1 a 2 ...a jo pa jQ+ i) ^Q) 

is a solution of the system (11) for j = j + 1. 

Proof is elementary and can be done by direct check. 

Lemma 1 given an efficient algorithm for construction of solutions of equations (11). 
In fact, solutions of these equations for j = 1 are well-known: they are Killing tensors 
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[21 EH> and a solution for arbitrary j may be obtained from a solution for j = 1 by 
successive application of the formula (26). If we manage to construct this way NJ 1 linearly 
independent solutions where N™ is given by the formula (25), then such solutions form a 
complete system in virtue of Theorem 1. 

Using the algorithm presented above we managed to obtain the general solution of 
equations (11) for m < 4 in the form 

pa^.-.aj _ giaj^^j pa 1 a2...a j - 2 ) _|_ ja x a 2 ...aj (27) 

where p a ^ a 2- a i-2 j s ^he general solution of equations (11) for j — > (j — 2) depending 
on iV™ 2 arbitrary parameters, and f a i a 2- a j j s a solution of equations (11) depending on 
NJ 1 — N™_ 2 arbitrary parameters. 

The first addend in the right-hand part of the formula (27) corresponds to such sym- 
metry operator (7) of order j that on the set of solutions of equation (1) can be reduced 
to a symmetry operator of order j — 2. Explicit expressions for f aia ^- a i corresponding to 
m < 4 are adduced below. 

1. m = 1. The corresponding tensor f a ^- a i can be reduced to a scalar not depending 
on the only variable. 

2. m — 2. Tensors f a ^ a 2- a i depend on two variables x\ and X2- The number of 
independent solutions, according to (25), is 

N = Nf - Nf_ 2 = 2j + 1. (28) 

Solutions are numbered by an integer number c satisfying the condition 

< c < j, (29) 

and include for c = one, and for each c > two arbitrary parameters giving independent 
components of a symmetric zero trace tensor \ a ^ a 2-- a j-c Q f ran k j _ c . The explicit form 
of the corresponding solution y" ia2 - a J j s given by the formula 

ja ia2 ... aj = £ ja ia2 ... aj + (1 _ £ )/>i a 2-%-i e aJ > Xfe , (30) 

where e aJb is the unit antisymmetric tensor, e = |[1 + (—1)°], 

ja\a 2 ...aj _ y(a\a 2 ...aj- c ^ y j rjfi 

/x=0 \i=j— c+1 

~.a 2k+l a 2k ) I (_ iyl L2-I ] (a^UfW* 



x 0- 




n^i h a a = v a ' - ' x 2 =x?+4 i=(-i^ (si) 

B < A, 



and symmetrization over the indices a\, a,2, ■ ■ ■ , a,- is implied. 
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3. m — 3. The tensor f a ^--- a o depends on three variables x — (x±,X2, x 3 ). The number 
of independent solutions is equal to 

TV = Nf - Nf_ 2 = i(j + l)(2j 2 + 4j + 3). (32) 

The solutions are numbered with pairs of integers c = (ci, C2) satisfying the conditions 

0<ci<2{|}, 0<c 2 <j-*{^±I}, ea = ^[l + (-l)°], (33) 

and include for each c the set 2ci + 1 of arbitrary parameters giving independent com- 
ponents of a symmetric zero trace tensor A aia2 "' ac i Q f ran k Cl _ Explicit forms of the 
corresponding solutions f£ ia2 ~- a 3 are given by the formula 

./:"" = ej™~ a * + (1 - e C2 )f b ^--- a ^e^ bc x c , (34) 

where e aJflc is the unit antisymmetric tensor, 

/A M +L M \ * /min({i},{i±i}-0 \ * 

= ^ ^Af"'^ [] X a *\ \\ g^k+A ^ 

p \i=A„+l / \ fc={ i (j 4 M+ L M )}+l / 

Here 

H = (^1,^2,^3,^4,^5), a; 2 = a; a x^ ab , 

5 M = 2 M2 + ^3 + /x 5 , ^ = 3 ~ ci - Z = (-1) C2+J ' +1 , 

L M = ci + /i 3 - 2/ii - /i 5 , F M = /i! + /x 4 , (36) 

and A Bm ' Am is an arbitrary symmetric zero trace tensor of rank + convoluted with 
vectors x^. 

Summation in (35) is to be done over all possible nonnegative values of /i satisfying the 
conditions 

0<^l<{|}, ^2 + ^3 + ^4={|}, < ^ 5 < ^[1 - (-1) C1 ]- (38) 

As well as in the formula (31) symmetrization over the indices ai, a 2 , ■ ■ ■ , cij in the 
right-hand part of (35) is implied. 

4. m — 4. The tensor j a ^- a i depends on four variables x = (x\, x 2 , x 3 , x 4 ). The 
number of independent solutions is equal to 

N = Nf-Nf_ 2 = i(j + 1)0' + 2)(2j + 3)(f + Sj + 4). (39) 
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Solutions are numbered with triples of integers c = (ci, c 2 , c 3 ) satisfying conditions (33) 
and (40): 



0<c 3 < j-2^^_±i K-2c 2 , (40) 



2 

and include for each c a set of iV c arbitrary parameters, where 
(c 2 + 2c 3 + l) 2 , ci = c 2 + 2c 3 



c \ 2(c 2 + 2c3 + l)(2 Cl -c 2 -2c3 + l), Cl ^c 2 + 2c 3 , 1 ; 

These parameters give independent components of an irreducible tensor 

yaia 2 —a Rl [a Rl+1 bi][a Rl+2 b2]~-[a Rl+R2 b R2 ] 

where R\ = c 2 + 2C3, i? 2 = ci — c 2 — 2c 3 (let us remind that an irreducible tensor of rank 
R\ + 2i? 2 has R\ symmetric indices and i? 2 symmetric pairs of antisymmetric indices, and 
convolution by any pair of indices and any triple of indices with completely antisymmetric 
tensor e^ pa vanishes). Explicit expressions for the respective solutions are given by the 
formula (42): 

j?aia 2 ...a,j _ ^^^/3 M ,(A M ,D C 



X 



fc={i(A M +D c +L M )}+l 




l [ g^k+l02 k ) I (42) 



where ji, x 2 , i\T M , are given by the formulae (36), (38) a,b = 1,2,3,4, = 

j - Cl - Bp - D c , D c = c 3 , I = (-l) n+1 , 

^_B M ,A M ,D C _ ^b 1 b2...b Bli aia2...a,A li [aA ll +idi]...[aA ll +D c dD c ] 
X X^ Xb 2 • • • Xfy B ^ Xd\ Xd.2 ' ' ' -Edjjc ) 

(43) 

symmetrization over the indices a±, a 2 , . . . , aj is implied in the right-hand side of (42). 

So, we have obtained the general solution of equations (11) in the space of dimension 
m < 4. One can verify by a direct check that the found solutions satisfy equations 
(11) and are linearly independent (it is not difficult to prove the latter considering i-fold 
convolutions of the found solutions with g kl ,0 < i < {§}). On the other side, these 
solutions form a complete system, as the number of arbitrary parameters they include is 
in compliance with the formula (25). 

Let us also mention that we can present the general solution of equations (11) also in 
the form 

j 

pa 1 a 2 ...a j _ yi 1 a 2 ...a l [aj+i&i] [ai +2 b 2 ]...[a j b j - l ] ^ ^ . . . Xb , (44) 

1=0 
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where \ a ^- a i[ a i+^\-\ a 3 b j-i\ j s a tensor, symmetric with respect to permutation of indices 
oi, . . . , dj and antisymmetric with respect to permutation of indices ai + j with bi, 1 < i < 
j — I, with vanishing convolution of this tensor over any three indices with e^ vpa . The 
latter means that a cyclic permutation with respect to any triple of indices (ak,ai +s ,b s ) 
gives zero, so the polynomial (44) admittedly satisfies equation (11). On the other side, 
the number of independent components of the tensor \ a ^--- a i[ a i+i b i\---[ a 3 b 3-i\ for < I < j 
is exactly N™ (25), so the formula (44) gives the general solution of equations (11). 
Decomposing tensors \ a ^--- a i[ a i+i b i\--\ a 3 b j-i\ ^ < I < j into irreducible ones (that is having 
vanishing convolutions with respect to any pair of indices), we come to formulae (26)-(43). 
We formulate the above results as the following theorem. 

Theorem 2. Equations (11) in a space of dimension m < 4 have N™ linearly independent 
solutions. These solutions are polynomials of x a of degree j and are given in explicit form 
by relations (26)-(43). 

The above theorem determines the explicit form of the Killing tensor of rank j in a 
space of dimension m < 4. 



5 Explicit form of symmetry operators Q n for n < 4 

The above results allow presenting in explicit form of symmetry operators of order n for 
equation (1) in m-dimensional space for arbitrary given n < oo and m < 4. For this 
purpose it is sufficient to look through all admissible values of c given by the formulae 
(29), (33), (40) and construct in accordance to the formulae (26)-(43) the corresponding 
expressions for Killing tensors of rank j p^ lCL2 --- a i (following (27), it is sufficient to restrict 
oneself with construction of f a ^- a j^ then substitute obtained expressions into (6), (7) 
and sum up over j from to n. 

In this Section we will realize this program for all n < 4 and m < 4, and write down 
in explicit form the corresponding symmetry operators. 

Let us calculate the number of linearly independent symmetry operators of order n. It 
is equal, according to (27) to the number of linearly independent solutions of the system 
(11) for j — n, n — 1 or (see (25)) 

! 2n 2 + 2mn + m(m -l)fn + m-2\fn + m-l\ 

N(n, m) = N"+ N2 = ; ^ , (45) 

v ' ; m m m(m-l) V rn-2 J\ m-2 J 1 v ; 

In particular, for m = 2, 3, 4 

iV(n,2) = (n + l) 2 , 

N(n, 3) = -(n + l)(n + 2)(n 2 + 3n + 3), 

N{n,A) = ^r(n+l)(n + 2) 2 (n + 3)(n 2 + 4n + 6). (46) 

f — 

Values of these numbers for n = 1,2,3,4 are adduced in Table 1. 



11 



Table 1. Number of symmetry operators of order n 
for equation (1) in m-dimensional space. 



m\n 


1 


2 


3 


4 


2 


4 


9 


16 


25 


3 


7 


26 


70 


155 


4 


11 


60 


225 


665 



Let us 
tions (11) 
m — 1 

n = 

m = 2 

n = 
n = 
n = 
n = 

n = 



write down explicitly the corresponding solutions F, 



1, 2, 3, 4, — X ln ; 



(n) 
m 



pai...a n Q f equa _ 



0. 
1, 

2, 
3, 



m = 3 

n = 
n = 
n = 



n = 



n 



^2 — A 20, 

^2 = A 21 + ^l^Xb', 

F a ia2 = g a ia2F (0) + A a ia2 + ^aigM)^ + X (g a ^ X 2 - X* 1 X* 2 ] 
F a ia2 a 3 = g {a ia2 pa,) + yaga, + ^i^gOs) 1 ^ 

+ A !o2)(f a2a3)x2 - a: a2 ^ aa) ) + \o,3)X {ai x a2 e a3)b x b ; 



(47) 
(48) 



/I 77101 020304 Aa\a 2 Tpa-idi) . ^01020304 , 1(010203 ,24) 6 
4, ^2 = 9 -^2 + A (0,0) + A fo,n £ 



W) 



+ \ { ^f(x a3 x a ^ - g a ^x 2 ) + \ {0A) (x^x a2 x a3 x a4) 
- 2x (ai a; a2 ^ a;ia4) + x 4 c/ (aia2 c/ a3a4) ); 



(49) 



0, Ff = A 3 ; 

1, F 3 a = A a + X b e abc x c ; 

n 7710102 r ,a\a 2 271(0) 10102 , \6(ai 

3 A (o,o) ~r A (o,r 

1 \b162 
A, 

771010203 Aa\a 2 T? a z) 1 \ 010203 

•3 ~~ f 3 A (0,0) 



\o,2) Xb i Xb 29 1 2 + A (1,0) X ^ 



X c -|- A/q 2)"^ ^^(0,2) X ^ X bi 

An 1 



- c a 2 )bc 
l (0,l) fc 



^(1,0) X 6i^' 



+ A ( 2,o)X ai x a2 

- A ( 2, )^ ia2 ^; 

771010203 „(aia 2 jp a 3) I \ 010203 I \*>(aio 2 03)60 1 \ 010203 2 _ 

r 3 — y r 3 -f a^ q ^ -f a^ q ^ t x c -f a^ q ^ x 

1 \ b i 6 2(ai o 2 a 3 ) , \&(aia 2 a 3 )6c 2 _ 91% 02 03)60 

(0,2) y dj bi' L b2 ' (0,3) x c x z,Aj-Q 2^ x t x c a 

_|_ )dbll>2 g(ai02 -«:()'/'• ,. , ,., \("l'i-> ,.« :! ) _ „«:>«:{) ,. j_ 



(50) 



, 1,6(0102 a 3 )6c 2 _ o\K«i 02 a 3 )6c 
>2 ~ (0,3) x c „l z,/\^q^^ „l c dj C dj(i 

.a 3 )dc , \(oio 2 a 3 ) _ \°( a i ^203) , a 2 a 3 )bc 

• h c Jj b\- L b 2 1 (1,0) (1,0)" " (1,1) 

A &d ^(0102^,03)60^ „. _ 1 \6 ^,(01^02^.03)60^ 



I \( a i ™0 2 03) 

-r Aj 2 jx x 
7701020304 



'(1,1) 

\i tl) g^ aia2 e a ' A)bc x c x d + Aj 2il) x (ai x a2 £ a3)fec 

\6 (0102^,03)60 2. 
2,0)^ /x (2,l)y c "°cX 1 

_ 20304 (01027711304) , \01020304 I \oi020304„2 , \6(ai0203 04)61 

3 -9^3+ A (0,0) + A (0,2) X + A (0,l) 8 

_ 91,6(010203 a 4 )„ _i_ \ fe ife(aia 2 0304)^ „ 1 16(010203 a 4 )6c 2 
(0,2) x -r a^q 2 ) y XfejX^ T A (o,3) fc x cX 
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x c 



n — 3, 



_ 9 x&d(oio 2 a 3 o 4 )6c , x 01020304 „4 , \d6i6 2 (ai a 3 a 2 a 4 )dc 

(0,3) Jj c-^d \ / '(o,4) (0,3) " *i-c* f 'Oi* t '02 

_ 4 \6(oi020 3 o 4 ) 2 ,4x6162(0102 o 3 a 4 ) r ^ ,93,6162(0102 a 3 o 4 ) r ^ 
(0,4) J/ 6 J/ ''(0,4) • l -bi J ^b 2 ' (0,4) " • b bi Jj 02 

- 4\^ 3{ai x a *g a ^x bl x b2 x b3 + A^ 3 *^ 

_ 3, K01 0203-04)60 , 3,(010203 04) -.2 _ 3,6(0102 0304) „ „2 

(1 1) ^ itj c it - / (i ~~ T~ (1 2^ (1 2) 

_ 9x6(0102 03 04) , qx 61 62(01 02 0304) 1 x(oio 2 o 3 „o 4 ) 

(1 "r (\ 2) y *^oi * j * y 02 "i - (2 0) 

_ 3,616263 -(0102 0304) _ x(oio 2 a 3 a 4 ) 2 , 3,6(01 a 2 a 3 a 4 )6c 

(1,2) x 6i^b 2 x 6 3 y y (2,0) " (2,1) c 

- \ 6<ai „020 3 -04)6c 2 , \(oia 2 a 3 04) 2 _ 9x6(01 a 2 a 3 a 4 ) 

(2 1)^* c<Xj \~ '\^2 2^ (2 2) ^ 

1 3,6162 (01 02 0304) _ 3,(0102 a 3 a 4 ) 4 , 9x6(01 a 2 0304) 2 

\^ ^2 2) 01 00 ^ i/ (2 2) y (2 2) y tJi'QtJU 

_|_ \( a l ™02 03 04) _ \6162 2 (0102-0304) _ \(oi 02 0304) 2 

^3 0) (2 2) ^1 ^2 y y (3 0) ^ 

- A^ 0) x 6l x (ai x a2 ^ a3a4) + A^ 0) x bl x 2 ^ (aia2 ^ a3a4) + A (4 ,o)X (ai x a2 x a3 x a4) 

- 2X m x^x a2 g a ^x 2 + A (4 ,o)X 4 ^ (aia2 ^ 3a4) ; (51) 



m = 4 

n = 0, F 4 {0) = A; 



», 1 I?o _ \a , x[odi] _ 

fl — ±, -T 4 — ^(0,0,0) "T A (0,0,l) Xrf i' 

_ paia 2 _ ^0102171(0) , xaio 2 , x(oi[a 2 di]) , x [oidi] [o2d 2 ] 

n — z, r 4 — g r 4 + A ( ,o,o) + A (o,o,i) Xd i A (o,o,2) ^di^efe 

, xaia 2 2 _ 93,6(01 02) 1 3,6162 aia 2 , 3,(01 o 2 ) 

' A (0,1,0) X ^ A (0,1,0) X X 6 "T A (0,l,0) X 6i Xb 2 y * A (1,0,0) X 

- Aj 1)0 ,o)^6^ aia2 + X mo) x a 'x^ - A (2 , ,o)^ ia2 x 2 ; (52) 

paia 2 a 3 _ (aia 2 rp a 3) , \aia 2 a 3 , X (0102 [o 3 di]) , x (01 [o2di][o 3 d 2 ]) 

x 4 —y r A ' (0,0,0) (0,0,1) Xd i A (0,0,2) X dl X d2 

x [aidi][a2d 2 ][o 3 d3] _9\6l(oio 2 o 3 )„ , X aia 2 a 3 2 

A (0,0,3) x dl x d2 x d3 ZA (o,l,0) x x 6j -r A (o,l,0) x 

3.6162(01 a 2 a 3 ) , x(oi02[o 3 rfi]) 2 _ 9 3,61 (01 [a 2 di] 03) 

(0,1,0) " X bl Xb 2 ~r ^(0,1,1) x x di Z,A (0,1,1) ^fei-^di 

, x6ib 2 ([oidi] a 2 o 3 ) 1 x(oio 2 a 3 ) _ 361(01 „a 2 a 3 ) , x(oi[a 2 di] 03) 

t A (0,1,1) y J-fei^6 2 x di -r A (i ; o,0) x A (l,0,0)y X 6j "T a (i i o,1) x x d x 

_ x6i([oidi] a2 o 3 ) 1 3,(01 02 03) _ 3,(01 0203) 2 , 3,([oidi] a2 03) 

A (1,0,1) i/ ■Hi-t'di ^ A (2,0,0) X X A (2,0,0)y X ^ A (2,0,1) X X Xrf i 

-AgS* 1 ^- >* 2 z dl ; (53) 

77^01020304 „(aia2 7710304) 1 \aia 2 a 3 a4 , x (aia 2 a 3 [a 4 di]) 

t A - g y t + A (0A0) + A (0 01) x dl 

x (ai02[o 3 di][o4d2]) , x (ai[a2di][a 3 d2][o4d3]) 

A (0,0,2) ^di^d2 T ^(0,0,3) •Ldi^d 2 x d :i 

, x ([aidi][o2d2][o3d3][a 4 d4]) , \ 01020304 2 _ 9 x 61(010203 04) 

r /\ Jj d\ Jj d 2 ' h dz' L d4, ~ (0,1,0) (0,1,0) 61 

, 3,6162(0102 0304)-. „ 1 x (oi020 3 [a 4 di]) 2 _ 0\6i(aiO2[o 3 di] a4 )„ „ 
(0,1,0) y x bi^b 2 T A (o,l,l) rf i (0,1,1) ^fei-^di 

1 x 6162(01 [a 2 di] 

) r y ^ rr. 1 \(oi02[o3di][o 4 d 2 ]) 2^ _ 
(0,1,1) y 61 "''62 -^di 1 A (o,l,2) ^di^d 2 

- 2X blMa2di][a;id2] x ai) x bl x dl x d2 + \ blb2i[aidl][a2d2] g a3ai) x bl x b2 x dl x d2 
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, 101020304 4 _ 4 \6l(oi020 3 04) 2 _|_4\ fe l 6 2(aia2 a 3 a 4 ) 

(0 2 0) (0 2 0) (0 2 0) • Xj b\ tlj b2 

2^(o,2,0) 9 3 4 ^ X bl X b2 ^ A (0,2,0) X 2 g 3 4 ^ Xj }1 Xb 2 Xb 3 
, n (aia 2 a 3 a 4 ) \brb2bab4 , \ 010203 ,^4) _ \&iaia2 0304)-, 

-|- y y A (o,2,o) x b 2 Xb 3 Xb 4 -t- ^(1,0,0) x A (i,o,o)y 

,(a 1 a 2 [a 3 di] a 4 ) _ \bi(ai[a 2 d 2 ] a 3 a 4 )„ „ , \ (ai[o2di][a 3 d 2 ] a 4 ) 
(1,0,1) rf i (1,0,1) " x fti x di ~r A (l 5 0,2) x ■ L d 1 - L d 2 

_ i6i(aidi][a 2 d2] o 3 a 4 )_ i _ _ ~ , \(«i«2«3 « 4 ) 2 _ ^(«i«2 a 3 » 4 ) 2 

(1,0,2) y xb 1 xd 1 xd 2 ~r ^(i^ o) x x A (i,i,o) y X bl X 

^ A (1,1,0) 4 ^ 4 ^°i ^ A (1,1,0) X 2 Q z 4 ^X bl X b2 ^(1,1,0)^^ 1 2 ^ 3 ^ x biXb 2 Xb:j 

I \( a i Q 2 a 3 a 4 _ \(ai«2 a 3 a 4 ) 2 , \(ai[a 2 di] a 3 „a 4 ) _ \(oi[o 2 di] a 3 a 4 ) 2 

A (2,0,0) X X A (2,0,0)y X A (2,0,l) X X Xd i (2,0,1) " X Xd l 

\([oidi][o2d 2 ] a 3 o 4 ) _ \([<ndi][a 2 d 2 ] a 3 a 4 ) 2 , >(aia 2 03 04) 2 

(2,0,2) -^di-^d 2 (2,0,2) " •^d\ dj d 2 ~ /v (2,l,0) 

_ r>\ (&i(ai a 2 a 3 a 4 ) , \fc162 (01 a 2 a 3 a 4 ) _ \(aia 2 a 3 a 4 ) 4 

(2,1,0) 6 i ~ l ~ (2,1,0) y x bi^b 2 /v (2,l,0)y 

, r>\bi(ai a 2 a 3 a 4 ) 2 _ \bib 2 (aia 2 ) a 3 a 4 ) 2 , \(ai o 2 a 3 a 4 ) 

~ 'V2 1 0) " (2 1 0)^ " ^ b\ iAj 02 > (3 0) 

— ~(oi~02 „a 3 a 4 ) _ \( a l a 2 a 3 a 4 2 , \bi „(aia 2 a 3 a 4 ) 2 

(3,o,o) x y Xfe i A (3,o,o) x y ^ A (3,o,o)y y 

_l_ \, s rr a \ rr a 2 rr a 'A ~<14 Q \ , ~(ai ~a 2 „a3a4) ~2 , \ (a 4 a 2 0304) 4 / cr/l \ 

+ A (4,0,0)£ X X X — ZA(4 i o,0)3 ; 2; y 'X + A(4 j o,0)^ <7 ^ • {^) 

Substituting (47)-(54) into (6), (7) and carrying differentiation operators to the right, 
we obtain explicit form of the corresponding symmetry operators. For n = 1 we have a 
complete set of symmetry operators of the following form: 

Q a 1 = P a = i-^, Qf = Jab = x a P b -x b P a . (55) 

We do not adduce explicit form of symmetry operators for n > 1 because of the 
corresponding formulae being extremely cumbersome (in fact, these expressions are given 
by relations (6), (7), (47)-(54)). 

First order symmetry operators Q\ adduced in (55), form a Lie algebra AP(p,q), 
satisfying the following commutation relations: 

[P a , P b \ = 0, [P a , J bc ] = i{g ab P c - g ac P b ), 

[Jab, Jed] = i{9acJbd + dbdJac ~ QacJbd ~ dbdJac)- (56) 

It is easy to notice using the representation (36) that symmetry operators of arbitrary 
order are polynomials of the operators (55). In other words, all symmetry operators of 
finite order of equation (1) belong to the enveloping algebra of the algebra AP(p, q). 



6 Explicit form of Killing tensor of arbitrary rank j 

Calculation of conformal Killing tensors of rank j (that is construction of the general 
solution of equation (14)) may be done similarly to what was presented above in Sections 
3-5. Construction of such solution simplifies utilization of the result formulated in the 
following lemma. 



14 



Lemma 2. Let p a ^--- a 3o be an arbitrary solution of the system (14) for j = jo, and F a 
be a solution of this system for j — 1. Then the function 

where [-] SL means the traceless part of the tensor in the square brackets (see (18)) is a 
solution of equations (14) for j — j + 1. 

Proof can be done by a direct check. 

We adduce below without proof the general solution of equations (14)for m < 4 and 
arbitrary j. 

By means of the reasoning similar to that in Section 3 we can show that in the two- 
dimensional space equations (14) are reduced to Cauchy-Riemann equations, and corre- 
sponding symmetry operators are determined up to arbitrary analytical functions deter- 
mining independent components of a symmetric traceless tensor F a ^ a -2-- a 3 (there are two 
such components for j ^ and one for j = (that is for the case when then tensor 
F aia2 - a i is reduced to a scalar). 

For m = 3 the number of independent solutions of equations (14) is equal to 

Nf = ±(j + l)(2j + l)(2j + 3). (58) 

Solutions are numbered by the pair of integers c = (ci, c 2 ) satisfying the conditions 

< ci < j, < c 2 < 2ci, (59) 

and for each c\ contain (2ci + 1) arbitrary parameters giving independent components 
of symmetric traceless tensor ,\ aia2 """ ac i f ran k; Cl Explicit form of the corresponding 
solutions is given by the formula 



p,aia2-..aj 

(cic 2 ) 



+ (i - e C2 )f\t:r 3 ' i£a3)bc ^ > ( c °) 



SL 



where 

/•Q,lQ,2---Q>j ^ " / C)\m 

*bib 2 ...bm(ai 

hac2) = 2^~^ \ m A (C 1C2 ) 

m=0 V 7 

X x a ci-m+i x a cl - m +2 _ _ . x aj) X bl X b2 . . . X bm X 2i{ ^ } ~ m) , (61) 

and the symbol [-] SL means the traceless part of the corresponding tensor; see (18) for 
m = 3. 

For m = 4 the number of independent solutions of equations (14) is equal to 

^ = ^(j + l) 2 (j + 2) 2 (2 J + 3). (62) 
The solutions are numbered by triples of integers c = (ci, c 2 , C3) satisfying the conditions 

0<ci<j, -c 1 <c 2 <c 1 , 0<c 3 < | Cl ~ |C2| |, (63) 
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and for each c contain N c arbitrary parameters where 

(ci + l) 2 , Ci = |c 2 |, 

2(|c 2 | + 2c 3 + l)(2ci - |c 2 | - 2c 3 + 1), ci ^ |c 2 |. 



(64) 



These parameters determine independent components of an irreducible tensor of rank 
R = i?i + 2^ 2 where 



Rt = |c 2 | + 2c 3 , i? 2 = ci - |c 2 | - 2c 3 , 
and explicit expressions for the corresponding solutions have the form 



(65) 



p(a 1 a 2 ...a j ) 



m+c 3 , 

i=0 v - 



m + c 3 w 2 



y ^6l&2---frm-i + c 3 (aia2-"a|c2|-m + i + C3 [ a | c 2 | + i- m+l + cg ■ • • ["cj -m + i-cg - | c 2 | -2c 3 ] 

x a .a cl _ m+i _ C3+ i a: ,a cl _ m+i _ C 3 + 2 . . . x ^) XhXb2 . . -x bm _ i+C3 



X X dl X d2 ■ ■ ■ £d ci _| C2 |_ 2c3 



(66) 



Here J\ 6l --- bm - i + c 3 a i--- a l c 2l- m + i + c 3[ a i c 2i+ i - m + 1 + c 3 dl l---[ ac i- m + i - c 3 a! ci-l c 2i- 2c 3] is an arbitrary irreducib- 
le tensor of rank Ri + 2i? 2 (i^ and R 2 are given in (65)), 



m 



-c 2 , c 2 < 0, 
0, c 2 >0, 



(67) 



( m ~j" C3 ) is a binomial coefficient, and the symbol [-] SL means the traceless part of the 
corresponding tensor; see (18), (19) for m — 4. Symmetrisation is implied over the 
indices a± , . . . , Oj in the righthand part (the sum over all possible permutations) . 

Thus, we have found the explicit form of the conformal tensor of rank j for m < 4. 
The formula (66) determines the general form such tensor for arbitrary m > 3, but at 
that the total number of independent solutions of equations (14) cannot be determined, 
in general, by the relation (62), but requires special calculation for each value of m. 

Let us point out that the general solution of equations (14) for m > 2 can be presented 
in the form 



paia,2...a,j 



j j-l-k 

^ ^ ^ ^ ybib2---b j -i-k-i{a\a2...a l+i [a l+i+1 d 1 ]...[ai +i+k d k ] 
Lk=0 i=0 



x b 



'j — l — k — i 



(68) 



where x b i b 2-bj-i-k-iai...a l+i [a l+i+1 d 1 ]...[a l+i+k d k ] j g a Censor S y mme tric with respect to permu- 
tation of the indices b\ , . . . , ai + j + k and antisymmetric with respect to permutation of the 
indices a;+j+/ with df,f= 1,2, ... ,k, with convolution of this tensor by any three indices 
with an absolutely antisymmetric vanishing. Decomposing such tensor into irreducible 
tensors we come to formulae (58)-(67) giving solutions of equations (14) for m = 3,4. 
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7 Examples of solutions and symmetry operators 
for n < 3 



Let us adduce an explicit form of the solutions obtained and corresponding symmetry 
operators for m < 4 and n < 3. Quantities of such solutions in accordance to (58) 
and (62) are adduced in Table 2. 

Table 2. Quantities of independent solutions of equations (14). 





1 


2 


3 


3 


10 


35 


84 


4 


15 


84 


300 



Quantities of the corresponding symmetry operators of order n can be obtained by 
summation of quantities of solutions from j — to j — n. We adduce the result in 
Table 3. 

Table 3. Quantities of order n symmetry operators of equation (1) 
with k = in m-dimensional space. 



rn\n 





1 


2 


3 


4 


3 


1 


11 


46 


130 


295 


4 


1 


16 


100 


400 


1225 



For m = 2 the number of solutions of equations (14) (and the number of the corre- 
sponding symmetry operators) is infinite as they are determined up to arbitrary functions. 

Explicit expressions for all independent solutions of equations (14) for m < 4 and n < 3 
are given by the following formulae (F^ = ftnia-aj)): 

m = 2 

3=0, = tp°(x 1 ,x 2 ); 

3 > <>• F 11 - 1 = (^- + ^)+<(^- + e;); 

F n - U = t(^-^)+^-Cr ( 69 ) 

Here ipj and £j are arbitrary analytical functions of two variables x±, x 2 , and other com- 
ponents of the tensor p a ^- a i are expressed through (69) using properties of of zero trace 
and symmetry. 
m = 3 

3=0, F^=X- 

3 = 1) ^(o,o) = \o,o)X a ; F(\^ = A° lj0) ; F ( " ;1) = e abc X b ^^x c ; 

F ( » 1)2) = 2X\ m x b x a - \ a {lt2) x 2 ; (70) 

3 = 2, ^ = A (0 ,o) (x a ^x a * - l -g a ^x^J ; 

rpa 1 a 2 _ \ai 02 1 \02 „ai _ ^ aia2 ^fe h 
^(1,0) ~~ A (1,0) X ~T A (1,0) X 3^ A (1,0) X ' 

7710102 _ / ai 02 1 T.02_ai N~6\c 
(1,1) — l X 5 6c + £ £ 60^^(1,1)1 
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j — 3, 



ciaia 2 _ / 01 \o 2 , a 2 \ «i \ 2 401 a 2 \b , aia 2 \ b 2 

^(1,2) ~~ y X A (l,2) + X A (1,2)) X 4X X A (l,2) X b T A (l,2) X b^ 



paia 2 
^(2,0) 



\ «l« 2 . 
A (2,0)' 



paia 2 ( a\bc\ba 2 , a 2 bc\6ai \ . 

^(2,1) — ^ fc A (2,l) fc A (2,l)J X o 

9 

pa x a 2 _ \aio 2 2 _ / 01 \a 2 b , a 2 \ aib \ , a\a 2 \bc . 

(2.2) — (2,2) V (2,2) T (2,2)/ b ~ ^(2,2) -H^cj 

paia 2 _ 9/ 01 u 1 ~o 2 ai \ „b /_ ai \a 2 fc , o 2 \aifc \ c 2. 

(2.3) ~~ V bk < X fc bfcJ A (2,3)- i x d l fc cfcA (2,3) cfc (2,3)/ 

paia 2 _ xoia 2 4 9/ ai \2c 1 T a 2 1,01c 2 , i 01 a 2 \cs 

(2.4) ~~ A (2,4) V X (2,4) ~r x A (2,4) x c x "T X A (2,4) x c^<i] 

1 



77x21 a 2 a 3 
^(0,0) 



paia 2 a 3 
^(1.0) 

rpaia 2 03 
^(1,1) 

771010:203 
^(2,0) 

771010:203 
^(2,1) 

rpaia 2 a 3 
pa\a 2 03 

*(2,4) 



A (0 ,o) ( x ai xV - -<? 



(0102^.03)^,2 



( a i „,a 2 „o 3 ) _ !/„( 

(i,o) x x gU/ 



= X 



(i,o) x2 + 

^x^^'A^d - ^g {aia2 E a:i \ c x b \ c ilA y, 

\(oio 2 03) _ _„(oio 2 \a 3 )6 
A (2,0) X A (2,0) Xfe > 

~(oi o 2 \0 3 )c _ 1 (oia 2 a 3 ) \cd b . 
^ c 6c (2,1) 6 bc / \(2,l) x x d> 

- 2x (ai x a2 A^ ) 2) x fe + -^ (aia2 x a3) Ag. 

^ 1 bc% 2 (2x 3 ^ A (2,3)"^rf A (2,3) X ) X ' 

x^Ag^x 4 - 4x a2 A^x b x 2 + 4x a2 x^Ag i4) x fc xO 



- -9 {aia2 { xaz) >HM) x kXiX 2 - \ a ^ b x b x 4 ); 



rpaia 2 a:j 
*(3,0) 
rpa\a 2 a 3 
^(3,1) 

771010:203 
^(3,2) 



^(3,3) 

■pa\a 2 a 3 
^(3,4) 



\ aia 2 a 3 . 

(3,0) ' 
£ (o3 bcA oio 2)6;Cc . 

» (010203) 2 _ 9^(03 \«i«2)d 1 ^_(oi02 \ a 3)bc . 
(3,2) C3.21 6 r " CS .21 6 c ' 



\3,2) 



(3,2) 



£^ *bc ^ A (3 5 3)^ x ^X 2 A^-g g^XftX c + 2 ^ A(3 3)XbX c Xrf 
\ (010203) 4 _ a (ai / \ a2i3)c 2 _ a 2 \ a 3)^c 



(3,4) 



- _^(«i^(2x a «)A^ 4) x d - Aj^x^x.x,; 

F oio 2 o 3 = £ (oi ( A g«)V - 4x a2 A^x.x 2 + 4x a2 x a3 )A^^x fc x / )x c 

™ia 2 a 3 _ \ (010203) fi _ c™(ai \0 2 03)c 4 
^(3,6) — A (3,6) X UX A (3,6) XcX 

4. 19 T (ai T «2 \ a 3) k Tl r r 2 — S r ( a i T a 2 r 03) \ 6cd 



m = 4 



7(0) 
(0,0,0) 



Ax a ; 



(71) 



(72) 
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^ 10) — \ % ccjj A x , 

F (a) _ \ [ad] 
r (1,0,0) ~~ A Xrf ' 

>) 
(1,1,0) 

-iaia2 

(0,0,0) ~ "(0,0,0) I x a, - -, 



41o) = A a ; (73) 
^(oTo) = Vo,o) (xV ~ ]g aia2 x 2 



pa 1 a 2 - 0\ b rf ( ^ a 2 _ \ n a\a 2 2 \ _ \ai a 2 2 

r (1,-1,0) ~~ ZA (1,-1,0) X & I x x ^i/ x j A (1,-1,0) X X 

_ \a2 „ai 2 | \ n a\a 2 \c 2. 

A (1,-1,0) X o y A (l,-l,0) x c x j 



F 



2* 

iaia 2 _ \[airfi] a 2 , \[a 2 di] ai . 
(1,0,0) A (l,0,0) x / \l,0,0) x 

^ rfl\0.2 \c _j_ \ai ™a 2 _j_ \a 2 „ai . 



paia 2 „aia 2 \c , \ai a 2 , \a 2 

r (1,1,0) ~~ 2 (1,1,0) c ' (1,1,0) tA (i,i,o) J ' ' 

paia 2 _9\bib 2 „,ai ™a 2 „. _ O \ fc i a i ™ a 2 rJ2 
(2,-2,0) ~~ A (2,-2,0) ^61^62 z,/v (2,-2,0) 6 l 

~~ 2 A (2,-2,0) X 1X Xf> i ^ A (2,-2,0) X 2^ 1 2 ^(2,-2,0) X X bi X b 2 'i 

F a ia2 _ x6l[aidi] ai , x&l^di] ai 

^ (2,-1,0) ~~ A (2,-l,0) x x &i x a!i ^ A (2,-l,0) x x 6i x di 

_ \ai[a 2 di] 2 _ \a 2 [aidi] 2 . 
(2,-1,0) d i (2,-1,0) d i' 

paia 2 _ \6oi a 2 , \ba 2 ai _ 9 \aia 2 2 _ l„aia 2 \bc . 
X (2,0,1) — A (2,0,1) X X& "l" A (2,0,1) X X& ZA (2,0,1) X 2 A (2,0,l) x & x c> 

paia 2 \[oirfl][o2d2] . 

1 (2,0,0) (2,0,0) x a!i^d 2 , 

j?aia.2 _ /\ai[a 2 d] \i2[aid]\ . 
X (2,1,0) ~~ V A (2,1,0) A (2,l,0) 'Xd-, 

r (2,2,0) ~~ A (2,2,0)' V V 

j = 3, F™ = A (0 ,o,o) (x^x^ - ; 

paia 2 a 3 _ \6 ( -(ai a 2 a 3 ) _ \ n (a\a 2 a 3 ) 2 ] _ \(<H 

r (1,-1,0) ~~ A (i,-i,o) I x x x 2 y / 6 (1,-1, 

+ ^ (aia2 ^ 3) A? l! _ 1 , 0) ^^ 2 + L^a^o/; 



^.02 ^,13)^,2 

Q-J it lAV 



6' 

/.'"J."2":i \ ,."2 ,.»:;) .. , _ 1,/ \ ,.2 . 

X (l,0,0) — A (1,0,0) X X Xd gi/ A (1,0,0) X Xd ' 

paia 2 a 3 _ \(ai a 2 03) _ 1 „(<Ha2 „o 3 ) \ 6 _ l„(aia 2 \ a 3) 2. 

r (1,1,0) ~~ A (1,1,0) X .^y X A (l,l,0) x 6 fi y A (1,1,0) X ' 

1 

— 9\ b ( ai T «2J3) T , T 2 1 n (aia 2 a 3 )\bc 2 , (aia 2 \ 13)^ 4 

(2,-2,0) 6 ^y (2,— 2,0)^^>' i 'C' i ' "T y (2,-2,0) b 

, x (aia 2 a 3 ) 4 , 1 (oia 2 \03)6 4. 
^ A (2,-2,0) X 3^ A (2,-2,0) XfeX ' 

F a ia2 a 3 _ xb([oid] a 2 a 3 ) \(ai[a 2 d] a3 ) 2 

^(2,-1,0) — A (2,-1,0) X X x b x d A (2,-1,0) X XrfX 



6' 

6" 

Z7iaia 2 a 3 _ xbife 2 / (ai a 2 a 3 ) _ ^ (bioj 03) 2 | 
(2,-2,0) — (2,-2,0) I 93 xxi XbjXfo 
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q9^ 1 2 ^(2^~l,0) X X b x d+ g</ 1 2 ^(2^-l,0) X X b X d', 

paia 2 a 3 _ \([aidi][a 2 d 2 ] a3 ) . 
(2,0,0) ~~ (2,0,0) X x di^d 2 , 

paia 2 a 3 _ ~a2~a 3 ) 1 (aia 2 a 3 )\6c 

^ (2,0,1) ~~ A (2,0,1) X X Xb 3^ X A (2,0,l) x 6^c 

4- I n K a 2 \ a 3)b 2 _ x(ai02 a 3 ) 2. 

+ A (2,0,1) X 6 X A (2,0,1) X X ' 

paia^ _ \(ai[<J2<i] a 3 ) _ 1 _(aia 2 \ [ a 3d])6 . 
^ (2,1,0) ~~ A (2,l,0) X Xd gf A (2,l,0) XfeXrf ' 

paia 2 a 3 _ \ (aia 2 a 3 ) _ ^ (aia 2 \ a 3)6 . 
(2,2,0) ~~ A (2,2,0) X 3^ A (2,2,0) Xb ' 

pa 1 a 2 a 3 _ \hb 2 b 3 | (01 02 03) , ^ „(aia 2 03) 2 

(3,-3,0) (3,-3,0) I -°oi J/ o 2 J '0 3 2" "i Jj o2 Jj o 3 

_ 1,(010303) 6 _ ox6i& 2 (ai a 2 o 3 ) 2 , o^6i(oio 2 a 3 ) 4 

(3 3 0) (3 3 0) 02 1 ^3 3 0) 01 

_ l fl (oio 2A o 3 )6ib 2 4. 
2 y A (3,-3,0) Xfc i Xfe 2 x J 

F a ia2 a 3 _ xbi&2([oid] a2 03) , * (aia 2 \ [o 3 o1)6ib 2 2 

-i ^3 _2 i 0) (3,-2,0) dj b\- L 'b 2 Jj d ~ gi/ (3,-2,0) dj o\- L 'b 2 Jj d 

_ 9X b( ai [a 2 d] o 3 ) 2 1 3, (aia 2 [a 3 d]) 4. 

__2 ^ ^b^d^ 1 (3 —2 0) ^d^ j 

pa ia2 a 3 _ \6i([aidi][o 2 d2] „a 3 )„ „ „ _ \ (aiNdi][a 3 d 2 ]) 2. 
^ (3,-1,0) ~~ A (3,-l,0) x x b^d 1 ^d 2 A (3,-l,0) -^di^d^ , 

pa ± a 2 a 3 _ \ [01 di] [a 2 d 2 ] [a 3 d 3 ] _ 
(3,0,0) ~~ (3,0,0) ^d 1 ^d 2 ^d 3 , 

T?a\a 2 a 3 /»6(oi[o 2 d] 03) _ I (aio 2 \ [«3d])&c \ _ \ (oia 2 [a 3 di]) 2. 

(3,0,1) ~~ ^ (3,0,1) X gi* (3,0,1) ^cj-^b^d A (3,0,l) Xd i X ' 

paia 2 a 3 _ f \ 6162(01 a 2 a 3 ) _ ^ (aia 2 a 3 ) \6i6 2 c 

r (3,-1,1) ~~ I A (3,l,l) X X 3^ A (3,-l,l) Xc 

1 l„(aia 2 \ 13)6162 „2\ _r ) \b{a 1 a 2 a 3 ) 2 , x(aia 2 o 3 ) 4. 

2 A (3,-1,1) X x bL^b 2 ZA (3,-1,1) X X 6 X (3,-1,1) ' 

paia 2 a 3 _ \ (ai[a 2 di][a 3 d 2 ] . 
(3,1,0) (3,1,0) ■ L d 1 ^d 2 i 

1 

Ciaia 2 a3 _ \ 6 ( a i a 2 a 3 ) _ _ _(aia 2 \a 3 )6c _ \(<Ha 2 a 3 ) 2. 

^(3,1,1) ~~ A (3,1,1) X Xfc ^iJ A (3,l,l) X 6 x c A (3,l,l) X ' 

pa x a 2 a 3 _ »(aia 2 [a 3 d]) _ 

(3,2,0) ~~ (3,2,0) X d' 

7-iaia 2 a 3 \aia 2 a 3 (7K\ 

r (3,3,0) — A (3,3,0) l'°J 

Here A(...), A" ^ A" 1 " 2 , A" 1 " 2 " 3 are arbitrary symmetric tensors with zero trace. 

Formulae (69)-(75) give explicit form of all linearly independent conformal killing ten- 
sors of rank j < 3 in spaces of dimension p + q = 2, 3, 4 (for p + q = 2 j is arbitrary). 
To obtain an explicit form of the corresponding operators it is sufficient to substitute 
(69)-(75)into (6), (7) that is write down j-multiple anticommutators of F aia2 - aj with 

d d d 
dx ai ' dx a2 ' ' ' ' ' dx a . ' 
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8 Killing tensors of rank j and order s 

Until now we considered solutions of equations (11), (14) that define Killing tensors (and 
conformal Killing tensors) of arbitrary rank j, but only of the first order. In this section 
we obtain explicit form of Killing tensors of rank j and of arbitrary order s. Such tensors 
are determined as general solutions of equations (16). 

The system of equations (16) is overdetermined including iV™ equations for iV™ un- 
known variables, where 



a™ I j + s + m- l\ - f j + m — 1 \ 



In the same way as it was done above in Section 3, we consider the set of differential 
consequences of the system under consideration that are obtained by £;-multiple differen- 
tiation of every term of the equation by jt~, jt—, • • • , jt— ■ This set is a system of linear 

OX a ^ OXa2 OXa^ 

homogeneous algebraic equations of the following form: 

pt(aia2--aj,aj + iaj + 2-..aj +s )bib2.-bk _ g (77) 

where the following derivatives are unknown variables: 

paia2...aj,aj + iaj +2 ...aj+sbib 2 ...b k _ Qaj+i QCtj+2 . . . QQj+s Qbi Qb 2 _ _ _ Qb k paia 2 ...aj (78) 

The quantities of unknown variables and of equations are equal to 

M k _ fj + m- 1\ fk + s + m- l\ k _ fj + s + m- 1\ fk + m- 1\ 
N »-{ m-1 ){ m-1 J' Nc ~{ m-1 ){ m-1 ) ' (?9) 

so conditions (24) are also fulfilled. 

It can be shown (see attachment) that the system (77) is non-degenerate, so it follows 
from (24) that 

piaia 2 ...aj,aj + iaj + 2...aj +a bib 2 ...bj _ g 

Whence we conclude that the Killing tensors of rank j and order s are polynomials of 
order j + s — 1. It follows from (77), (79) that such polynomial contains arbitrary 
parameters, where 

s-l / . . -, \ / . . - \ s-l 

m 



i—n \ / \ / h—n 

(80) 



i=0 x ' x 7 fc=0 

s fj + m — 1\ fj + s + m— 1 
m \ m — 1 J \ m — 1 

Here the first sum gives the number of independent solutions that have order by x 
smaller than s. Such solutions can be written in the form 

F a ia2 ... aj = \£g-g x bi x i» ... x b % t<s, (81) 
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where A^ 2 "^ are numeric parameters with no limitations set by equations (77) (certainly 
the symmetry with respect to permutations of indices a\ <-> a' x , b^ <-> b' A, A' = 1, 2, . . . , j, 
= 1,2,..., i. 

In the case s = 1 the formula (80) is reduced to (25). In particular, for m = 2,3,4 we 
obtain from (80) that 



2 
1 

12 



n % = ~^<3 + 1)0' + 2)0 + s + 1)0 + s + 2), 



4 = ^s(j + l)(j + 2)(j +3){j + s + l)(j + s + 2)(j + s + 3), (82) 

Thus we have determined the number of linearly independent Killing tensors of rank 
j and order s. To compute these tensors in explicit form we will use the following two 
lemmas. 

Lemma 3. Let p^ ia2 --- a ^ fa Killing tensors of rank j and order s, and ip be a functions 
satisfying the equation 

d> 1 d u ip = 0, n,v=l,2...m. (83) 

Then the function 

Paia2...aj _ ^pa-Laz-aj j-g^-j 

is a Killing tensor of rank j and order s. 

Proof is reduced to direct check of the lemma statement that is (s + l)-multiple 
differentiation of (84) by d 0- ^ 1 , d aj+1 , . . . , (9 aj + s+1 and subsequent symmetrization of the 
obtained expression by a\, 02, Oj+s+i using relations (16), (83). 

Lemma 4. Let p a ^--- a ] he a Killing tensors of rank j and order s. Then the convolution 

paia,2~-aj-\ _ paia,2- ■■ a j % (85) 

is a Killing tensor of rank j + 1 and order s + 1 . 
Proof is similar. 

The adduced Lemmas provide an effective algorithm for construction of Killing tensors 
of order s from Killing tensors of order 1 found above in Section 4. The only difficulty 
in application of this algorithm is the need to sort out all linearly independent solutions 
of the system (16) (whose number is determined by the formulae (80), (82), as, generally 
speaking, there are more solutions of the form (84) than we need). 

The general solution of equations (16) is determined in the following theorem. 

Theorem 3. Equations (16) in the space of dimension m < 4 have n™ linearly indepen- 
dent solutions where n™ is given by the formula (82). These solutions have the form 

jpaia2—a.j _ (aj-iaj rp a i a 2---aj-2) , ia 1 a 2 ...aj 
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1 

+ ^2x a i +1 X a i+ 2 ■ ■■X a i+ d - l) , Ej = -[1 + (-l) j ], (86) 
d=l 

where F aia2 --- a 3+<i-i are Killing tensors of rank j + d — 1 and of order 1 whose explicit form 

(s) 



is given by Theorem 2, p, aia2 - a ^ 2 are Killing tensors of rank j — 2 and of order s; 



■2.-1 

( 3 - — 1\ 



x ^+2^+3 . . . ^■_ 2 a,_ 1A [a j ),c] ;rc) ^ 

^[aj,c] ^ s an ar l)it rar y antisymmetric tensor of rank 2. 

Proof. By virtue of Lemmas 3, 4 the function f^° 2 -°j given by formula (86) is a Killing 
tensor of rank j and order s; the first term — by definition, the second — in accordance to 
Lemma 3 (being the product of of a Killing tensor of order 1 and tp = X^x^), the third - 
in accordance to Lemma 4 (each convolution with x M lowers the rank and increases the 
order of a Killing tensor, and we make the first convolution with the first-order tensors 
described above). 

It is to some extent more difficult to make sure that formula (86) gives all linearly 
independent Killing tensors of order s. Proof of linear independence of all terms of the 
formula (86) is reduced to comparison of terms having the same order by x ai using different 
convolutions by one, two etc. pairs of indices. Calculation of then umber of independent 
solutions given by formula (86) can be done easily by sorting through independent solu- 
tions for first-order tensors F aia2 '" a i+ d - 1 entering the last term (such tensors are described 
in Theorem 2, but it is necessary to restrict consideration with solutions corresponding 
to c 2 > d — 1, as others give zero input into convolutions of (86)), and adding the number 
of solutions of the form (87). The result obtained is in compliance with the formula 

where n™ s is the total number of solutions given by the formulae (82), n™_ 2s is the number 
of solutions of the form g( a i- ia i p^ a2 --- a J- 2 ) j s a j so gj ven j n (82), iV is the total number 
of solutions under the summation sign and of solutions of the form (87). We omit the 
corresponding cumbersome calculations. 

Formula (86) determines recurrent relations for calculation of explicit form of a Killing 
tensor of rank j and order s from a known tensor of order s and rank j — 2. Such 
calculations may be easily checked starting from known killing tensors of order 1 and 
arbitrary rank, see Theorem 2. 

Let us adduce as an example explicit expressions for Killing vectors of order s < 3 in 
three-dimensional space received from general relations (86): 



s = 


1, 




= A a -I 


- e abc r] b x c , 




s = 


2, 


rpa 


= F h 


+ \ ab x b + \x a + e abc r] bd x c x d + £ V - 


x a i b x b 


s = 


3, 


rpa 
*(3) 




+ X abc x b x c + \ b x b x a + e abc r] bdl x c x d x l 
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Here e abc is a unit antisymmetric tensor, and other Greek letters designate arbitrary 
symmetric tensors with zero trace. 



9 Conformal Killing tensors of rank j and order s 

Let us briefly discuss equations (17) describing conformal Killing tensors of rank j and 
order s, and adduce without proof solutions of these equations for arbitrary j, s and 
m < 4. 

A constructive way for finding solutions of equations (17) is shown by the following 
statement that may be proven by direct check. 

Lemma 5. Let ^r® 102 -°j ^ e a conformal Killing tensor of rank j and order s, and ip be 
an arbitrary function satisfying the equation 

d»d v y = g^X, A = const. (88) 

Then the function 

is a conformal Killing tensor of rank j and order s + 1 . 

We can show that quantities of linearly independent solutions of equations (17) for 
arbitrary j, s and m = 3, 4 are given by formulae 

m = 3, Nf s = S -{2j + l)(2j + 2s + l)(2j + s + 1), 



6 

s 

12' 



m = 4, Nf s = ^(j + l) 2 (j + s + l) 2 (2j +2 + s). (90) 



Using Lemma 5 we managed to construct iV™ linearly independent Killing tensors of 
rank j and order s (giving full system of solutions of equations (17)) in the following form: 

Fr ,.. B , = £ f Fr ...a j{x2) ^ tfA . (91) 

i=l V d=0 J 

Here p^ ia2 - a ^ are conformal tensors of rank j and order 1 given by formulae (60), (61) 
or (66) (the index "z" distinguishes independent solutions of (91) with various degrees 
of x 2 ), fili 2 d " a3 are tensors of rank j whose explicit form is adduced below. 

In the case m = 2 the number of conformal Killing tensors of order 1 appears to be 
infinite, see (69). The same formulae (69) give a general form of a conformal Killing tensor 
of arbitrary order s. 

In the case m = 3 functions f^ 2 d ' aj are characterized by an additional integer c, 
< c < 2j, (92) 
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and are determined up to arbitrary symmetric zero-trace tensor \ a ^-- a R Q f rank R = 
2(c + i) — 1. Explicit form of these functions is given by by formula (93): 



/■a\a2...dj 
id 



£ Ji d T'" aj + (! - £ C )f b{aia2 - a i- 1 e aj)bc 



.1 ( 



SL 



(93) 



where 



aia,2...aj 
idc 



n=0 



X X 



Q>2 • • •Q'j — n 



l j—n+l 



X 



n 



Clj-n+2 . . . T a j) 



X 3 X^Xbz • • • Xb d+n X 2 



2({§}-n) 



(94) 



with the symbol [•] designating the zero trace part of the corresponding tensor, see 
(18), (19) for m — 3, and the index d is introduced for numeration of linearly independent 
solutions of (91) with different degrees of x 2 . 

In the case m = 4 the functions f°^ 2 ' aj are characterized by a pair of additional 
indices c = (ci, c 2 ) 



-J < ci < j, < c 2 < 



J - Ci 



(95) 



and is determined up to arbitrary irreducible tensor ,\ oia2 "" OH i [afli+i&i]"[ a fli+fl 2 b «2] f rank 
i?i + 2R 2 where 



-Ri = |ci| + 2c 2 + «, R 2 — j - |ci| - 2c 2 . 
Explicit form of these functions is given by the formula (97): 



(96) 



a\a,2---Q>j 
id 



n+C2 / _|_ \ 

( — 1)^ / 71 2 J ^ rc 2 ^(i^6i62---fen+ C 2-d-i(aia2---a| cl |-„ + d+i+ C2 +irfi]- 



■[ a j-n+i+d-C2^j- |ci | — 2c 2 ] j. a j - n+i+d- c 2 + 1 . . ■ X " 3 Xh ' ' ' Xh Xrl ' ' 'Xrl 

"1 "n+C2— d — i ^1 < ^j — \ci\—2c2 



SL 



, (97) 



where 



n 



-ci, ci < 0, 
0, ci > 0, 



(98) 



and symmetrization is implied in the right-hand part by indices oi, a 2 , . . . , %. 

The formulae (91)-(98) give in explicit form all linearly independent conformal Killing 
tensors of rank j and order s in space of dimension m = p+q < 4. In particular, conformal 
vectors of order s < 3 in three-dimensional space, in accordance to (91), (94) have the 
following form: 



s = 1, Ffo = \ a {1) + e^^Xc + ? {1) x 2 - 2x a e {1) x b + fix a ; 
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s — 2, — F^ + F^x 2 + A^2)^6 + £ abc i](2) x c x d + C(2) a ' 2a '6 — 2x Q ^| ) 2)XfeX c ; 

S = 3, = -F(2) + £ 4 -F(i) + x2 (^l3) x b + £ a C V(3) x c x d + C^J) 3 ^ 2 — 2x a £;| ) 3)X;,X c ) 

_i_ \ abc T, _i_ r abc bdk , afcc 2 _ 9 a^bcd 

(3) & c "i - ^ '/(3) -Lc-Ld-Lk "T C( 3 )X£,X C X ZX ^(3) X^XcXd- 

Here the unit antisymmetric tensor, and other Greek letters designate arbitrary 

symmetric traceless tensors F^, F^ and are first-order Killing vectors (generally 
speaking, different). 

10 Conclusion 

Let us sum up. We have defined the notion of Killing tensor of rank j and order s and 
of conformal Killing tensor of rank j and order s. These tensors are defined as general 
solutions of equations (16) or (17) that in the case s = 1 coincide with generally accepted 
equations for Killing tensors and conformal Killing tensors, see e.g. |llj . 

We limit ourselves with investigation of equations (16) and (17) in flat de Sitter space, 
and generalization of these equations for for the case of spaces with non-zero curvature 
requires replacement of d ai for covariant derivatives. 

Equations (16) and (17) are natural generalizations of the Killing equations and 
arise in description of higher-order symmetry operators. In the present paper we show 
relation of these equations (for first-order tensors) with higher-order symmetry operators 
of Klein-Gordon-Fock equation, see Sections 2, 5. Equations for Killing tensors and 
conformal Killing tensors of order s < 1 arise in problems of description of symmetry 
operators of order s for systems of partial differential equations - - in particular, for 
the Maxwell equations [13]. We have found in explicit form all non-equivalent Killing 
tensors of rank j and order s in the space of dimension p + q for arbitrary j and s and 
p+q < 4. Limitation by dimension of space is based, on the one hand, on practical reasons 
(the absolute majority of equations of mathematical physics being the field of research 
interests of the authors, have dimension m < 4 with respect to independent variables), 
and, on the other side, on difficulties that had not been overcome to the moment in proof 
of non- degeneracy of systems of algebraic equations for coefficients of Killing tensors in 
spaces of arbitrary dimension, see attachment. At that formulae (42), (66), (86) giving 
solutions of equations (16), (17) for p + q = 4, probably give the general solution of these 
equations for arbitrary p + q > 4. 

The found general solutions for Killing tensors and conformal Killing tensors of arbi- 
trary order and rank may find quite large use in description of symmetry operators of 
systems of partial differential equations. In this paper using these solutions we found full 
set of symmetry operators of arbitrary finite order for Klein-Gordon-Fock equations with 
zero and non-zero mass. 
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A Non-degeneracy of systems 

of equations for coefficients of Killing tensors 

We will adduce proof of Theorem 1 stating non-degeneracy of system of linear algebraic 
equations (22). As it is cumbersome we adduce it in abridged form. 

The main difficulty of the analysis of system (22) is the need to do it for arbitrary value 
of j, that is for a system of arbitrary fixed dimension given by formula (20). 

Let us consider equations (22) for m — 4, at that equations for m < 4 will be included 
into the analysis as particular cases. Indices a±, a 2 , ■ ■ ■ , %+i and bi, b 2 , • • • , bk with k < j 
independently take values from 1 to 4. At that, as it is easy to notice, the system (22) 
splits at non-linked subsystems M(si, s 2 , S3, s 4 ), where si (I — 1, 2, 3, 4) gives the number 
of indices having the value /. It is obvious that 

si + s 2 + s 3 + s 4 = j + 1 + k, (A.l) 

so < si < j + k + 1. 

System (22) is non-degeneracy iff all its subsystems M(si, s 2 , s 3 , s 4 ) are non- degeneracy. 
Without loss of generality, for arbitrary subsystem M(si, s 2 , S3, s 4 ) we can put 

si < s 2 < s 3 < s 4 , (A.2) 

other cases can be reduced to (A.2) by renumeration of variables. 

Let us prove non-degeneracy of an arbitrary subsystem M(si, s 2 , S3, s 4 ). 

We designate by the symbol n\ (I = 1,2,3,4) the number of indices of unknown variable 
pa 1 a 2 ...a j ,a j+1 b 1 b 2 ...b k p resen ^ on ^ e [ e ft Q f comma an d equal to I, and by the symbol 

mi — the number of indices after the comma that are equal to I. Obviously, the following 
should be satisfied, 

iti. c + n c = s c , n\ + n 2 + n 3 + n 4 = j, rri\ + m 2 + m 3 + m 4 = k + 1, (A. 3) 

so out of eight numbers rii and rrii only three will be linearly independent (see (A.l)). Let 
us choose the following numbers as independent: rii, n 2 and n 3 , then the triple (rii, n 2 , n 3 ) 
will completely determine a vector j? a i a 2---a-j,a j+1 b 1 b 2 ...b k f rom ^g subsystem M(s±, s 2 , S3, s 4 ). 
Using for such vector the designation F(ni, n 2 , n 3 ) and considering relations (A.1)-(A.3), 
we can write any equation (22) from the subsystem M(si, s 2 , S3, s 4 ) in one of the following 
forms: 



(n 2 + 1)F(0, n 2 , n 3 ) + n 3 F(0, n 2 + 1, n 3 - 1) + (j - n 2 - n 3 )F(0, n 2 + 1, n 3 ) = 0, 
max{0, Si + s 2 — k — 1} < n 2 < s 2 — 1, 



(n 3 + 1)F(0, 0, n 3 ) + (j - n 3 )F(0, 0, n 3 + 1) = 0, 
max{s 3 — fc — 1, —1} < n 3 < s 3 — 1; 



(A.4) 



max{0, Si + s 2 + s 3 - - 1 - n 2 } < n 3 < min{s 3 , j - n 2 }; 

(ni + l)F(ni, n 2 , n 3 ) + n 2 F(ni + 1, n 2 - 1, n 3 ) + rz 3 F(ni + 1, n 2 , n 3 - 1) 



(A.5) 



+ (j - ni - n 2 - rz 3 )F(ni + 1, n 2 , n 3 ) = 0, 
max{0, Si — fc — 1} < rii < Si — 1, max{0, s x + s 2 — fc — 1 — r^} < n 2 < s 2 
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max{0, si + s 2 + s 3 — k — 1 — n\ — n 2 } <n 3 < min{s3, j — n\ — n 2 }. (A.6) 
When s 1 > (the case Si = is considered below) there are three possibilities: 



1. si + s 2 < k + 1, 

2. si + s 2 > k + 1, si < k+ 1, 

3. si>Jfc+l. 



(A.7) 



Let us consider these possibilities one by one. 

In the case 1 the system under consideration is given by the formulae (A.4)-(A.6). 
Let us present the vector F(ni, n 2 , n 3 ) in the form of a column whose components are 
numbered by the index 



F(n 1} n 2 , n 3 ) = (F(n 1} n 2 , h 3 ), F(ni, n 2 , n 3 + 1) . . . F(ni, n 2 , n 3 )) T , 



(A.8) 



where n 3 and ri 3 are minimal and maximal values of n 3 , and each vector F(n 1 , n 2 ,h 3 + k), 
in its turn, will be regarded as a column whose components are numbered by the index n 2 , 
< n 2 < s 2 . Then it is possible to write equations (A.4)-(A.6) in the matrix form: 



AF = 0, 



(A.9) 



where 



A = 



B, 



( B 

Ei B\ 

2E 2 B 2 

\ siE B ai J 

( A 

En A+i 

2Ei 2 Di +2 



\ 



\ 



(A.10) 



s 2 Ei S2 Di +S2 J 



Here and Eif (I, k — 1, 2, . . . , si, f — 1, 2, . . . , s 2 ) are unit matrices whose number of 
rows coincides with the number of rows of the adjacent matrices on the right (the number 
of columns of B k (D{) coincides with the number of rows of B k+i (D i+1 )), Di being the 
matrices whose explicit form is determined below. Namely, with s x + s 2 + s 3 < k + 1: 



f .1-1 + 1 



D, 



\ 



3-1 
2 j-l-1 



s 3 j - I + 1 - s 3 J 



(A.H) 
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with s\ + s 2 + s 3 > k + 1 

/ a\ 3 — ai — I + 1 



ai + 1 



j — ay — l 



\ 



\ 



a 2 3 -a 2 -l + l J 



(A.12) 



where ai = max{0,Si + s 2 + s 3 — k — /}, a 2 = min{s3, j — / + 1}. 

At that in the cases a± = or (and) S3 > j — I + 1 in the matrix (A.12) the first or 
(and) last column should be crossed out. 

Note A.l. It can be shown that there are always will be less of matrices Di of the 
form (A. 11) (or (A.12) ai ^ j - a x - I + 1 ^ 0) than of the matrices (A.12) with 
ai = j — a 2 — I + 1 = 0. 

Our task is to prove that all rows of the matrix A (A. 10) are linearly independent. 
Writing this matrix in the equivalent form 



A' 



( E 1 B x 

2F 2 B 2 



\ B 



\ 



siE sl B Sl 



J 



and subjecting A' to the transformation A' 



rank, with V and W being reversible matrices of the form 
/ E 1 



V 



A" = VA'W that does not change the 
\ 



E 2 



E 



Sl 



\ -B B B 1 ... (-iy^J p B B 1 B 2 

I Ei —Bi ^B\B 2 —^BiB 2 B 3 . . . 
E 2 —B 2 



B x , En 



w 



\\E 2 B Z 



-V Sl T^iy.BiB 2 . 



B Sl -i \ 



^ S1 - 1 JlhyB 2 B 3 ...B. 



E, 



si 



-B 



Sl-l 



si-1 



E. 



si 



we get 



/ Ei 



A" 



Eo 



E 



Sl 



\ 



jpB B 1 B 2 ■ ■ ■ B S1 J 



(A.13) 
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and proof of linear independence of the rows of the matrix (A. 10) is reduced to proof of 
linear independence of the rows of the matrix 



B Sl = B Bi ■ ■ ■ B Sl . 



(A. 14) 



Lemma A.l. The matrix (A. 14) can be split into blocks Bfj including the same number 
of rows as the matrix D Sl+ i^i, 1 < I < s 2 +i, 1 < / < $2+1- These blocks have the following 
form: 



B lf 



0, -s 2 < I - f < or si + 2 < I- f < s 2 , 

Sl + *) p/r/A-iA • • • D si+f _ u < 1 - f < Sl , 



i-f 

si + 1 

Si + 1 



(A.15) 



P l -'E, l-f = Sl + l. 



Here D1-1, Di, . . . are matrices from (A. 10), the general form of which is given by the for- 
mulae (A. 11), (A. 12), and E k are zero and unity matrices of corresponding dimensions, 



pk 



k\ 



k\(m-k)\' m~ (k-m)\ " 

Proof can be done by induction, by successive investigation of the products B B 1 , 
B BiB 2 , . . .. 

Lemma A. 2. The matrix (A. 14) by finite number of elementary transformations can be 
reduced to the form 



( 



B 



D D 1 . . . D S1+S2 \ 



D sl -iD sl . . . D S2+ i 
D Sl D Sl+ i ■ ■ ■ D S2 



E 



S2~S\ 



Eo 



\E, 



(A.16) 



/ 



Proof. We can describe elementary transformations mentioned in the lemma in the 
following way: 

1. Let us go from the matrix (A. 14) to the matrix Tl d (with d successively taking values 
1, 2, . . . , s 2 — Si), using the following algorithm: 

1) Represent the matrix (A. 14) in the block form (A.15) and operate with "rows" with 
the number I, including blocks B^ , B^, ■ ■ ■ , Bfl 2+1 . 

2) Multiply (t + l)-th "row" of the matrix B Sl by the matrix D t _i, and the "row" of 
the matrix B S1 with the number t — by the number Si — 1 + 2 and subtract the latter from 
the former. Write the obtained result instead of the "row" with the number t. Perform 
this operation successively with all "rows" for t — 1, 2, . . . , si. 

3) Multiply (si + 2)-th "row" of the matrix B Sl by D Sl . . . D S2 . . . D si+d _i and subtract 
from the obtained result the "row" with the number s\ + 1 multiplied by d. Write the 
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obtained result instead of the "row" with the number si, leaving the remaining "rows" 
unchanged. 

4) As a result of the described transformations B Sl — > U' d , where IT^ is the matrix 
having in the first column the only non- vanishing element (II^) Sl+21 . It is possible to get 
by means of elementary transformations that all elements of the Si + 2-th "rows" would 
also vanish (except (n^) Sl+2 i)- 

5) Let us put the s\ + 2-th "row" to the lowest position. As a result we get the matrix 



n. 



1 - S2 of the following form: 
/ 



n 



S2-S1 



E, 



Ei 



M \ 



J 



(A.17) 



where M is a matrix that can be split into blocks of the following form: 



0. 



si<l-f<0, 



s 2 + l 
l-f 



Pl:{D l „ 1 D r --D f+S2+1 , 0<l-f<s 



(f,l) = l,2,...,s 



i+i- 



2. Let us simplify the matrix M by using successively the adduced algorithm for 
q = 0, 1, . . . , si — 1, and for each value of q — for k = 1, 2, . . . , si — q. 

1) Let us multiply the (k + l)-th "row" of the matrix M by the matrix Dk-\, and the 
k-th column — by the number s\ — k + 2, and subtract the former from the latter, leaving 
other "rows" unchanged. 

2) Perform this operation successively for all k — 1, 2, . . . , s\ — q, simplifying the matri- 
ces obtained at each step by means of elementary transformations vanishing all elements 
of "rows" except one that was the only non-vanishing in its column. 

3) Perform operations 1), 2) successively for all q = 0, 1, . . . , si — 1. 
As a result we come to the matrix (A. 16). The lemma is proved. 

Lemma A. 3. Let the matrix D = ||<i a ft||, a — 1,2, ... ,s, b = 1,2 ... ,r, s < r have the 
rang s, with all minors D being positive, and the matrix B have the form 



( b 



B 



ii 
b 2 \ 



\ 



b 2 2 
b 32 



'r-lr-2 



V 



J r—lr—l 
Kr-1 ) 



(A.18) 



where bkk > and bt+ik > 0, k = 1, 2, . . . , r — 1. Then the matrix C = DB also has the 
rank s, and all minors C are positive. 

Proof is reduced to direct utilization of the Binet-Cauchy formula [TS| representing 
spinors of the matrix product DB via sum of the products of minors of the matrix D by 
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minors of the matrix B. As a result each minor of the matrix C can be represented as 
the sum of positive values. The Lemma is proved. 

By virtue of the Lemma A. 2 proof of linear independence of the matrix A (A. 10) is 
reduced to proof of linear independence of rows of the matrix D d , 



D Sl -dD Sl -d+i • ■ ■ D s 



d = 0,l,2,...,si, 



(A.19) 



D Sl -d+i, ■ ■ ■ are matrices of the form (A. 11) or (A. 12). By virtue of Note A.l 



V d 

where D s 

the number of rows of each matrix (A.19) does not exceed the number of its columns. 
Considering successively the products D sl _ d D Sl ^ d+ i, D Sl _ d D Sl _ d+ iD Sl _ d+2 , . . . and using 
each time either the Silvester inequality 16 or Lemma A. 3, it is not difficult to show 
that the rank of the matrix (A.19) coincides with the number of its rows, and, whence, 
all rows of the matrix A (A. 10) are linearly independent. 

We have proved non-degeneracy of of the system (A.4)-(A.6) for the case 1 from (A. 7). 
In the case 2 when s% + s 2 > k + 1, the system (A.4)-(A.6) is reduced to equations (A. 5), 
(A. 6) that can also be written in matrix form (A. 9) where A is given by (A. 10), but the 
blocks B\ have a new form. Namely, with s 2 < k + 1 



I {s 1 + s 2 -k- 1)Eq 



V 



[Si + s 2 



Do 

-k-l 



l)£i Dj 



s 2 E, 



k— s\+l 



\ 



(A.20) 



0, 1, . . . , si + S2 — k — 1, and 

^ Dl+k-Si~S 2 

E 1 D 



B, 



l+k— si— S2+1 



V 



(A.21) 



■-■2 



if Si + s% — k < I < Si. If s 2 > k, then all matrices Bi, < I < si, are given by the 
formula (A.20). 

Explicit form for the matrices Di is given by relations (A.22)-(A.23): 
for s 2 < k 

I CL\ j — I — d\ \ 

a 1+1 



D, 



j - I - ai - 1 



V 



(A.22) 



a 2 j - I -a 2 J 

where a± = max{0, Si + s 2 + s 3 — k — I — 1}, a 2 = min{s 3 , j — I}, I — 0, 1, 
/di j — l + k + 1 — si — s 2 



. k — di, 



D, 



l + k + 1 

cii 



si 
1 



3 -l + k 



\ 



s l — s 2 ~ a l 



a 2 



j - I + k - si - s 2 - a 2 J 
(A.23) 
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where a± = max{0, S3 — I}, a 2 = min{s 3 , j — l+k+1 — s± — s 2 }, I = k—si + 1, k — si+2, . . . , k. 

If ai = then first columns in (A. 22) and (A. 23) should be crossed out, if s 3 > j — I, 
then the last column in (A. 22) should be crossed out, and with s 3 > j — I + k + 1 — si — s 2 
it is necessary to cross out the last column in (A. 23). 

In the case s 2 > k the matrices Di are given by the formula (A. 22) for all I. 

Our task is to prove linear independence of rows of the matrix A determined by re- 
lations (A. 10), (A.20)-(A.23). Transforming this matrix to the form (A. 13) we reduce 
this problem again to investigation of the matrix (A. 14) that in our case can be split into 
blocks of the form 



B Sl 



6, — s 2 <l — f<Si + s 2 — k or k — s 2 + 2 < I — f < s 2 , 

\l + s 1 + s 2 -k-f) P ^ + s,-k- f A-i A • • • A^ 2 -i-/, 



Sl 



Si 

1 



Si + 1 



s 2 + k < I - f < k - s 2 , 

E, l-f = k + l 



(A.24) 



S2, 



where 1 < / < k + 1 - s 2 , 1 < / < s 2 + 1. 

Further proof is done in full analogy with the proof for the case 1. 

Let us consider now the third case from (A. 7). The corresponding system of equations 
(A.4)-(A.6) is reduced to equations (A. 6). Writing these equations in the matrix form 
(A. 9) we come to the corresponding matrix A of the following form: 



A = 



where 



( ( Sl -k)E B \ 

(s 1 -k + l)E 1 B 1 

\ siE k B k J 



{s 2 - R)E Q D 

(s 2 -R+l)E 1 Di 



(A.25) 



R = 0,l,...,k, 



D, 



\ s 2 E R D R ) 

f (s 3 - I) j + k + 1 - si - s 2 - s 3 

s 3 - I + 1 j + k - si - s 2 - s 3 



V 



S3 j + k + 1 - si - s 2 - s 3 - IJ 



I 



0,1,...,*;. 

All rows of the matrix A (A.25) are evidently linearly independent. 

Thus we had proved that the system of equations (A.4)-(A.6) is non-degenerate in 
all cases listed in the formulae (A. 7). Thus the system (72) in the case m = 4 is non- 
degenerate. 
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Considering only such systems of equations (A.4)-(A.6) that correspond to Si — we 
get a full set of non-linked subsystems of the system (72) for m — 3, and in the case 
si = s 2 = we come to full set of non-linked subsystems of the system (72) for m = 2. 
Consequently non- degeneracy of the system (72) for m = 2 and m = 3 follows from the 
adduced proof particular case. 

Similarly (but with involvement of somewhat more cumbersome calculations) it is pos- 
sible to prove non- degeneracy of the systems of linear algebraic equations (74) for coeffi- 
cients of Killing tensors of rank j and order s. 
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